
for data = [20, 50, 100]

  clf;
  points = strcat('points_', num2str(data));
  allpts = strcat('gp_data_', num2str(data));
  D = load(allpts);
  dim = size(D);
  n_class = max(D(:,dim(2)));
  P = load(points);
  [C,ia,ib] = intersect(P(:,1:2),D(:,1:2),'rows');

  for c = 1:n_class
    S = D(ib, :);
    vals = S(find(S(:,3)==c),:);
    s = strcat('@', num2str(c));
    plot(vals(:,1), vals(:,2), s);
    hold on;
  end
  pause;

end



% colours are indexed rgbmcw